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Progress in computing inclusive B decay spectra 

Einan Gardi and Jeppe R. Andersen 

Cavendish Laboratory, University of Cambridge, 
J J Thomson Avenue, Cambridge, CB3 OHE, UK 

We review the progress in the QCD calculation of inclusive decay spectra. It was recently shown that the inher- 
ent infrared finiteness of inclusive spectra extends beyond the level of logarithms. Dressed Gluon Exponentiation 
makes practical use of this property by computing the Sudakov exponent as a Borel sum. Based on renormalon 
analysis, infrared sensitivity in the exponent is reflected in power corrections that are inversely proportional to 
the third power of the mass. Therefore, the parametric enhancement of non-perturbative corrections near the 
phase-space boundary is effective only in a small region. Consequently, the on-shell decay spectrum provides a 
good approximation to the meson decay. In particular, it facilitates a precise determination of \V u b\ from present 
measurements of inclusive charmless semileptonic widths without involving a non-perturbative "shape function" . 



One of the important challenges for QCD pre- 
sented by the B factories is the calculation of 
inclusive partial decay widths within specific re- 
gions of phase space. The obvious example is the 
inclusive measurement of charmless semileptonic 
decay, B — ► X U W, used for the determination 
of | V u b | • Owing to the overwhelming charm back- 
ground, this measurement is strictly restricted to 
the region of small hadronic mass, Mx < 1.7 
GeV, where charmed final states are kinemat- 
ically excluded; specific cuts within this region 
are adopted depending on the experimental tech- 
niques applied. Consequently, extracting \ V u b\ re- 
lies on quantitative theoretical understanding of 
the spectrum. 

The common lore has been that estimates of 
the partial widths for the relevant cuts strongly 
depend on the non-perturbative structure of the 
meson. The spectrum in the small— -Mx region 
was obtained [1-3] through the convolution of 
a computable hard coefficient function with a 
leading-twist momentum distribution function of 
the b quark in the meson, the "shape func- 
tion" [4, 5]. Being non-perturbative, the lat- 
ter was not computed, but rather parametrized 
and fitted to the measured B — ► X s j spec- 
trum, in analogy with deep inelastic (DfS) 
structure function phenomenology. The "shape 
function" has been the dominant source of 
uncertainty in determining \V u b\. 



Recently the situation has changed, ft was 
shown that the inherent infrared finiteness of in- 
clusive spectra extends beyond the level of log- 
arithms [6]: the leading renormalon ambiguity 
in the Sudakov exponent cancels against that 
of the pole mass while the next-to-leading one 
is absent. Upon using Dressed Gluon Expo- 
nentiation (DGE) [7-12] the B — > X s j and 
B — ► X u lv decay spectra can be well approx- 
imated by the corresponding on-shell decay spec- 
tra [6,13-17]. Non-perturbative effects associated 
with the meson structure enter in this framework 
as power corrections, and have just a small impact 
on the experimentally-relevant partial widths. 
The purpose of this talk is to explain why this 
is so. 

Inclusive B decays such as B — > X s j and 
B — > X u lD are dominated by jet-like momen- 
tum configurations: the hadronic system X has 
a small mass, Mx <C Mb, while its energy in 
the B rest frame is large. Inclusive decays can 
be analyzed within perturbation theory owing to 
the fact that the heavy quark carries most of the 
meson momentum and it is therefore close to its 
mass shell. By making the perturbative, on-shell 
approximation one neglects non-perturbative ef- 
fects that are suppressed by inverse powers of the 
heavy-quark mass. In the region of interest, of 
small Mx, these power corrections are paramet- 
rically enhanced. This will be a central issue in 
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what follows. 

Given the typical jet-like momentum config- 
uration it is convenient to compute the spec- 
trum in terms of lightcone coordinates. The 
hadronic lightcone coordinates are defined by: 
P ± — E x T \Px\- The jet is characterized 1 by 
large P~ , of 0(Mb) and small P + , which is of 
the order of the mass difference between the me- 
son and the quark 2 , A = Mb— mi,. Within the on- 
shell approximation, we define the partonic light- 
cone coordinates: pf = P^ 



the p~j distribution is trivial 



- A. At Born level, 
5(p~j). Beyond this 
order, the small p^ region is characterized by soft 
and collinear gluon emission that broadens the p^ 
distribution. This is why Sudakov resummation 
is essential in computing the spectrum. 

One obvious limitation of perturbation theory 
is the fact that the perturbative spectrum, at 
any order, has support only for p^ > 0, i.e. for 
P + > A, while the physical spectrum fills up the 
whole region down to P + = 0. As we shall see 
below, this inherent limitation is removed in the 
DGE approach [15]. 

In perturbation theory, the large hierarchy be- 
tween and pj in the peak region is reflected 
in large Sudakov logarithms. It is convenient to 
choose the resummation variable as the exponent 
of the rapidity, r = P^/pJ and write the triple 
differential width in b — > X u lT> as [17] 



1 dT(X,r, xi) 
Tr, dXdrdxi 



V(\xi)d(r)+R{X,r, Xl ), (1) 



where the additional kinamatic variables arc A = 
pj /nib and xi = 2Ei/mb, the Born-level width 
is T = G F \V u b\ 2 ml / (192ir 3 ) , the virtual correc- 
tions are 



V(X,xi) = w (X,xi) 



wi(A, xi) + 



(2) 



and the real-emission contribution, which starts 
at 0{a s ), is further split as R(X,r,xi) = 



1 In B > X s j the photon is real so P = Mb, while 

in B > X U W, where the lepton pair has a mass q 2 = 

(M B - P + )(M B - P~), the distribution in P~ is rather 
broad but still peaks near P~ ~ Mb, see Fig. l4l 
2 Here and below mi, stands for the quark pole mass. The 
renormalon ambiguity [18,19] is dealt with explicitly using 
the Principal Value prescription. 



R Biae .(X,r,xi) + R res XX,r,x t ). The regular piece 
-Rreg. (X, r, Xi), similarly to V(X, x{), is known [17, 
20] to 0(a s ) only, while 3 
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R sine XX,r,x t ) = -w (X,xi)[ h-- 

r 4 r I „ tt 



which contains all the non-integrable terms at 
r — ► 0, regularized as ()» distributions, is known 
in full to 0(a 2 s ) — see Eq. (3.41) in Ref. [17] - 
and in part 4 , at higher orders. This higher-order 
information is contained in the resummation for- 
mula: 



dT N (X, xi) 



1 



(l-xi)/X 



_ — / dr(l ., 
dXdxi To Jo dXdrdxi 

= V(X, xi) x Sud(m 6 A, N) + AR N (X, x t ). (3) 

In the small p^ region, corresponding to large N, 
the first term dominates: it sum up to all orders 
perturbative corrections that diverge as powers 
of IniV or are finite in the N — ► oo limit. The 
remainder, 



r) 



at_i dT(X, r, xi) 



AR N (X,xi) 



(l-x,)/A 

drO.-r) tf - 1 R ns Xr,^,xi), 



is of 0(1/N) and can be included at fixed order. 

The resummation formula J3Jl is a manifesta- 
tion of the infrared and collinear safety of the 
on-shell decay spectrum: owing to the cancel- 
lation of infrared singularities between real and 
virtual corrections, the moments have finite ex- 
pansion coefficients to any order in perturbation 
theory. As usual, this cancelation leaves behind 
large Sudakov logarithms. These arise from two 
distinct subprocesses [6, 21-23]: the formation 
of a jet of invariant mass squared of Q(pfpj) 
and radiation off the nearly on-shell b-quark, 
characterized by transverse momenta of 0(p~j'). 
The two corresponding functions, the jet func- 
tion Jm [12,24-26] and the quark distribution in 
an on-shell heavy quark Sn [14, 27, 28] obey the 



3 For the definition of the ()» distribution and the details of 
the NLO result in these variables see Sec. 3.1 in Ref. [17]. 
4 The exponent in Eq. below is known to NNLL accu- 
racy. Nevertheless, in order to determine i? 3ing . (A, r, a;;) 
to this accuracy at O(a^) and beyond one would need to 
compute the NNLO virtual coefficient W2(X, x{). 



3 



following evolution equations in moment space: 

(4) 

where f = 1 — r and where and S are 
Sudakov anomalous dimensions that are now 
known in full to NNLO [14,15,28,29]. The Su- 
dakov factor, summing up the logarithmically- 
enhanced terms in decay spectra to all orders, 
is: Sud(Amb, N) — Jjv(Ami; /XF)5jv(Am&; //f), 
where the factorization-scale dependence cancels 
exactly in the product. 

Note the fundamental difference with DIS 
structure functions where there is just one 
source [12] of Sudakov logarithms, namely the jet 
function Jn(Q, (J-f), and the factorization-scale 
dependence cancels against the non-perturbative 
quark distribution function q^dip), e.g. ~ 
Jn(Q', /J-f)qn(^f)- Contrary to Sn, <?jv cannot 
be defined in perturbation theory owing to the 
collinear singularity of an incoming light quark. 

Because of the cancellation of logarithmic in- 
frared singularities, the evolution kernels defined 
by the r.h.s. in {2} are finite to any order in per- 
turbative theory. However, these kernels conceal 
infrared sensitivity at the power level [6-17, 30- 
32], which only becomes explicit once running- 
coupling effects are resummed to all orders [7-12]. 
A systematic way to quantify this infrared sen- 
sitivity is to regularize the ensuing divergence 
of the perturbative expansion by Borel summa- 
tion [33-35]. To this end one writes the scheme- 
invariant Borel representation [36] of the anoma- 
lous dimensions [6-17]: 



J{a s {^)) = ^-j du[^\ T(u)Bj(u), (5) 
Po Jo \M / 



2 \ M 



00 Jo 



du 



T(u)B s {u), 



where /3q 



\N f and T{u) = 
{uS) uS e- uS /T(l +'uS) with S = PJffi. Using © 
in one can explicitly perform the r integra- 
tion and identify potential power like ambigui- 



ties arising from the r — ► limit as Borel singu- 
larities. The solution of the evolution equations, 
formulated as a Borel sum, is [15]: 



Sud(Arribi N) 
x B s (u 



exp 



flo Jo 
T(-2u)T(N) 



' u 



A 2 



X 2 m 2 



-Bj{u) 



1 

2a 



2u) 

r(N - u) 



(6) 



In the jet-function part one finds potential renor- 
malon ambiguities at positive integer values of u, 
while in the soft-function part at any integer 
and half integer u. Being anomalous dimensions, 
Bj(u) and B$(u) are not expected to have any 
renormalon singularities of their own, however, 
unless these functions vanish at these locations 
there will be power-like ambiguities in the evo- 
lution kernels in (@J. These ambiguities scale 
as powers of A 2 N/(Anib) 2 and AN/(Xmb), cor- 
responding to the jet-mass and the soft scale, re- 
spectively, so they are parametrically-enhanced 
at large N. 

The discussion of renormalon singularities can 
be made concrete by focusing on the gauge- 
invariant set of radiative corrections correspond- 
ing to the largc-/?o limit 5 . The anomalous di- 
mensions are then obtained as analytic functions 
in the Borel plane [12]: 



es 



1 



1 



sin itu 



Bj{u)\ x&xse0o = — [j_ u , I _ u/2J nu 
Bs(u)\ laIsef3o =ei"(l-u). (7) 

Based on these expressions one can deduce which 
power ambiguities indeed appear in Eq. ©. On 
the soft scale one finds ambiguities correspond- 
ing to u = h and u > | and on the jet mass 
scale, u = 1 and u = 2. Of course, these ambigu- 
ities should all cancel once non-perturbative cor- 
rections are systematically included. Conversely, 
in absence of a non-perturbative calculation, the 
ambiguities provide a good hint on the functional 



5 Results in this limit are obtained [34] by first considering 
the large— Nf limit, in which a gluon is dressed by any 
number of fermion-loop insertions, and then making the 
formal substitution Ne > — 6/3q. 
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form, and quite possibly even the magnitude, of 
the non-perturbative power corrections. 

Although the dominance of non-perturbative 
corrections that are probed by renormalons can- 
not be established from first principles in QCD, 
this conjecture has led to successful power- 
correction phenomenology in a variety of appli- 
cations [34,37-39]. The DGE approach was ap- 
plied to several infrared and collinear safe observ- 
able, notably event-shape distributions [7, 8] and 
heavy-quark fragmentation [10, 11] where precise 
data from LEP was used to test the renormalon 
dominance assumption. The observed hadroniza- 
tion effects near the phase-space boundary were 
found to be in agreement with the pattern of 
power corrections predicted by renormalons. On 
these grounds one expects that the renormalon 
structure of the Sudakov exponent provides a 
good indication on non-perturbative corrections 
also in B decay spectra. Of course, here the lead- 
ing corrections are associated with the initial- 
state B meson, rather than final-state hadroniza- 
tion, and the renormalon dominance assumption 
must be again confronted with data. 

In computing the Sudakov factor © we assume 
that the pattern of renormalon singularities of the 
large-/?o limit holds in the full theory. Nonethe- 
less, the perturbative expansions of these func- 
tions at u — and likewise the residues in © get 
modified by O(l//?o) contributions. In order to 
match {7J onto the perturbative expansions of J 
and S [14], we write the following ansatz [15,17]: 



x e ClU+c ^" 2 x Wj(u) (8) 



Bj(u)=Bj(u)\ laIsef3o 

B s (u)=B s (u)\ laigcPo x e Cl "+ c ^ 2 x W s (u), 

where ci = (l - tt 2 /4) C a /(3/3 ) is univer- 
t is 

sal [27], c 2 are determined based on the 
known NNLO expansions of these functions and 
Wj/s(u) parametrize yet-unknown 0(u 3 ) cor- 
rections. By fixing Bs(u) and Bj(u) this way 
and choosing the Principal Value prescription 
for the renormalons, the Sudakov factor © is 
uniquely determined. 

Now, the resummed spectrum is obtained by 
an inverse Mellin transformation of Eq. @ , 

1 dT(\,r,xi) f dN 



To dXdrdxi 



2ni 



{1 -ry N ^^l, (9) 
dXdxi 



NLO+LL 

0.4 NLO+NLL 

NLO+NNLL 

DGE 

0.3 - 




P + 



Figure 1. The P + spectrum in B — > X u lv in the 
on-shell approximation at a representative point 
in phase space: P~ = 4.79 GeV; E t = 1.71 GeV. 
The spectrum is computed with fixed logarithmic 
accuracy (LL, NLL, NNLL) and by DGE. All re- 
sults are matched to NLO [17]. 



where the integration contour C runs parallel to 
the imaginary axis, to the right of the singularities 
of the integrand. Finally, the triple differential 
distribution in physical, hadronic variables is: 



1 dT(P+, P~, Ei 
T dP+dP-dE, 
1 dT(X,r,xi) 



(10) 



To dXdrdxi 



P+-A 



The change of variables in (|10fl has a subtle but 
absolutely crucial role in obtaining the correct 
distribution within the on-shell approximation. 
To understand it observe that this transformation 
involves A = Mb — rnj, that is ambiguous since 
the pole mass mi, has a u — ^ renormalon [18, 19], 
and recall that the Sudakov exponent in also 



has a 



renormalon ambiguity. These ambi- 



guities have the same source namely the ambigu- 
ous definition of an on-shell state and they cancel 
out exactly [6] in Eq. i|10|l . This cancellation was 
checked explicitly in the large-/?o limit and it is 
understood to be general. In the course of the 
calculation we deal with it by using the Principal 
Value Borel sum regularization in both Eq. © 
and in the calculation of the pole mass [15] in A. 
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Figure 2. The average energy in B — ► X s j, 
(i? 7 ), in the on-shell approximation, as a func- 
tion of the minimum photon energy cut Eq, as 
calculated by DGE in [15] varying m" s within its 
error range. The result is compared with data 
from Belle [42]. Inner and total error bars show 
systematic and statistical plus systematic errors 
(added in quadrature), respectively. 

The extent to which the perturbative on-shell 
spectrum is constrained depends primarily on 
having the soft anomalous dimension function 
Bs{u) entering Eq. JfjJ) sufficiently well under con- 
trol. Eq. (JSJ incorporates the QCD perturba- 
tive expansion of Bg(u) around the origin up 
to NNLO. However, owing to the soft scales 
probed as N gets large, which are 0(Xmb/N), 
the Borel sum has some sensitivity to Bs(u) 
away from the origin. This sensitivity depends 
on what is assumed about Ws(u) in Eq. ©. 
Here one can make use [15] of additional infor- 
mation, namely the value of the leading renor- 
malon residue of the pole mass 6 , which fixes 
Bg(u — |) owing to the exact cancellation of 
ambiguities explained above. We further assume 
that Ws(u) has no Borel singularities and thus 
the vanishing of Bs\i aise ^ (u — 1) carries over 
to the full theory. It was shown [15, 17] that 
so long as Bs(u) does not get particularly large 
beyond this region, e.g. near the u = | renor- 

6 The pole-mass u = i renormalon residue has been accu- 
rately determined, see Refs. [15,40,41]. 



Figure 3. The energy variance in B — > X s j, 
({(Ej) — -E 7 ) 2 ), as a function of the minimum 
photon energy cut Eq, as calculated by DGE 
in [15] varying a s within its error range, com- 
pared with data from Belle [42]. 

malon position, the Principal Values Borel sum 
in Eq. © is well under control. Moreover, the 
resulting spectrum in hadronic variables H10[l has 
approximately the correct physical support prop- 
erties, i.e. it smoothly extends into the non- 
perturbative region P + < A and vanishes for 
P + ^ 0. This highly non-trivial property of 
the Borel-resummed on-shell spectrum suggests 
that non-perturbative corrections in this frame- 
work are not large. 

Having obtained this result through resum- 
mation, it is interesting to return to the con- 
ventional Sudakov resummation framework, in 
which Eq. JHl is expanded and computed to a 
given logarithmic accuracy, and see why it fails. 
The comparison between the results is shown in 
Fig. n Obviously, there is a qualitative dif- 
ference: the fixed-logarithmic-accuracy curves 
are sharply peaked reaching a Landau singular- 
ity near P + = 0.5 GeV where they become com- 
plex; in contrast, the DGE result, which is not af- 
fected by Landau singularities, smoothly extends 
to the non-perturbative regime P + < A ~ 0.4 
GeV. Aside from the Landau singularity issue, 
one observes that results of increasing logarith- 
mic accuracy represent a badly divergent series, 
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which is dominated by the leading u = i renor- 
malon. The divergence sets in early (see Tables 
1 and 2 in Ref. [16]) and therefore increasing the 
formal accuracy in this naive approach may re- 
sult in worse approximations. In contrast, in the 
DGE approach the renormalon is regularized and 
the associated ambiguity cancels. 

Eventually, the dependence of the on-shell de- 
cay spectrum on the adopted Principal Value pre- 
scription concerns renormalons at u = k/2, with 
all integer k > 3. These give rise to ambigu- 
ities scaling as powers of (AN/(\m b )) K . The 
corresponding power corrections can be summed 
up into a non-perturbative "shape function" , 
T (AV/(Am&)), which multiplies the Sudakov fac- 
tor. This function has a clear field-theoretic in- 
terpretation 7 as the moment space ratio between 
the quark distribution in the meson and that in 
an on-shcll heavy quark [6]. A priori, power cor- 
rections on the soft scale could affect the entire 
peak region. However, based on the renormalon 
ambiguities we actually find the situation is quite 
different: the absence of the first two powers and 
the inherent suppression of higher powers that 
is dictated by the structure of Sudakov expo- 
nent (Eqs. (3.50) and (3.51) in Ref. [17]) sug- 
gest that the shape function should mainly affect 
high moments N > Xnib/A and therefore be im- 
portant only in the close vicinity of the endpoint, 
P + — ► 0. It amounts to small corrections else- 
where. 

Given the properties of the resummed spec- 
trum and the relatively minor role of non- 
perturbative corrections in this approach, one can 
directly use the perturbative on-shell result as an 
approximation to the meson decay spectrum; in 
this approximation T (AiV/(Am&)) = 1. This was 
suggested in Ref. [15] where predictions for the 
B — > A s 7 spectrum were obtained. In addition 
the first few central moments with a varying lower 
cut on the photon energy were computed. Soon 
after, first results for these moments were pub- 
lished by the BaBar collaboration, which agree 

7 Note that this is quite different from what became the 
standard terminology in the B-decay community, where 
the term "shape function" is used as synonymous to the 
quark distribution in the meson, denned with some ultra- 
violet cutoff. 



T 0-35 L u 




P' [GeV] 



Figure 4. The P~ spectrum in B — ► X U W as 
calculated by DGE [17], after integration over P + 
and Ei in four different situations: no cuts (full), 
Ei > 1 GeV cut only (dotdashes), E t > 1 GeV 
with P+P- < M% = (1.7 GeV) 2 (dashes), and 
Ei > 1 GeV with P+ < P+ ax = 0.66 GeV (dots). 



well with the predictions — see Fig. 4 in [16]. 
A similar comparison with Belle results [42] was 
done later on and is shown in Figs. |2] and |3J 

The good agreement [16] of the on-shell cal- 
culation with the measured B — ► X s j decay 
spectrum provided a strong incentive to apply 
this approach to B — > X U W. This was done in 
Ref. [17] where \V u b\ was extracted for the first 
time directly based on resummed perturbation 
theory, without relying on any parametrization of 
the spectral shape. As explained above, the re- 
summation applies to the fully differential width. 
This means that the measurement of any partial 
branching fraction can be readily translated 8 into 
a measurement of \V u b\. 

Fig. 0] presents the computed event-fraction as 
a function of P~ . It shows the distribution for 
two cuts that have been used by Belle [43] to dis- 
criminate charm: an upper cut on the hadronic 

8 The DGE calculation has been implemented numer- 
ically in C++ facilitating phase— space integration with 
a variety of cuts. The program is available at 
http : //www . hep . phy . cam . ac . uk/~andersen/BDK/B2U/] 
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Figure 5. The event fraction in B — > X u lv 
within the range P + P~ < M x and Ei > 1 GeV, 
as calculated by DGE, plotted as a function of the 
cut value Mx- Different assignments for the un- 
known parameter C (see Sec. 3.2 in [17]), which 
is proportional to B$(u — |), are used. 

invariant mass at Mx =1.7 GeV, and an upper 
cut on P + , at -P+ ax = 0.66 GeV. In both cases 
a lower cut on the lepton energy E\ > 1 GeV is 
applied as well. Belle data in these measurements 
are: 

AB(P + P- < (1.7 GeV) 2 ) = 1.24 • 10 -3 (±13.4%) 
AB{P + < 0.66 GeV) = 1.10 • 10 -3 (±17.2%), 

where the numbers in the brackets represent the 
total experimental error. The DGE results for the 
corresponding event fractions are: 

Rc Ut (P + P~ < (1.7 GeV) 2 ) = 0.615 (±9.6%) 

R C ut{P + < 0.66 GeV) = 0.535 (±15.2%). 

Using 

AB(B — > X u lv restricted phase space) 

= r B r tot (B — » XJD) x i? cut , (11) 

and T tot (B — -4 XJD) = \V ub \ 2 x 66.5 ± 4p S -\ 
we obtain [17] 

\V ub \ = (4.35 ± 0.28 [cxp] ± 0.26 [th] ) • 10- 3 
\V ub \ = (4.39 ± 0.36 [cxp] ± 0.40 [th] ) • 10- 3 , 



respectively. The largest source of uncertainty, in 
both the calculation of the total width and in that 
of the event fraction, is the value of the short- 
distance mass m" s . Apart from a precise mass, 
to further improve this determination it would be 
necessary to have a complete NNLO perturbative 
calculation to the fully differential width. Impor- 
tantly, we find that the sensitivity to the details 
of the quark distribution function is small. An es- 
timate of this source of uncertainty was obtained 
in Ref. [17] by changing Ws(u) in Eq. Fig. [5] 
shows the result as a function of the Mx cut. 
Evidently, the effect is small for experimentally 
relevant cuts. 

Needless to say, the potential of the DGE ap- 
proach is far from being exhausted by this pertur- 
vative determination of \V u b\. In order to quan- 
tify the leading power corrections that constitute 
J- (KN / {\mb)) one would obviously need careful 
comparison with the measured spectrum. Both 
B — ► X s j and B — ► X u lv decays can be used 
for this purpose. 

To conclude, we have shown that by properly 
resumming the perturbative expansion, the on- 
shell approximation provides reliable predictions 
for inclusive B decay spectra, facilitating a precise 
determination of |V^b|. The calculation of the Su- 
dakov exponent as a Borel sum guarantees renor- 
malization and factorization scale invariance, and 
opens the way for incorporating valuable informa- 
tion on the large-order behavior of the series. At 
the end of the day, the main advantage of this 
approach stems from the possibility to use the 
inherent infrared safety of decay spectra, which 
extends beyond logarithmic accuracy. 
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